Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

金属银熔体快冷过程的计算机模拟

利用分子动力学模拟技术研究了由周期性边界条件控制的500个银(Ag)原子构成的金属Ag熔体快速冷凝过程.冷却速率为8×1013K/s.模拟在FS(Finnis-Sinclair)相互作用势的基础上,通过双体分布函数、键对分析技术、键取向序等多种方法,对液银快速冷凝过程的微观结构转变特性作了分析,给出了连续快速冷凝过程中液银原子间依靠相互作用力形成的独特的微观结构图像.并考察了冷却过程中体系能量和元胞体积随温度的变化.模拟结果表明在快速冷凝过程中液Ag没有形成bcc结构的倾向.     

锥形光纤中光波传输特性的蒙特卡罗模拟

为了分析锥形光纤中光子流的传输特性和空间分布特性，采用了蒙特卡罗算法对其进行了模拟仿真，得到了大量光子经过锥形光纤后的光子分布图.模拟实验结果表明，通过锥形光纤的大量光子集中在探测器中心，到达探测器某一位置的光子数随着该位置到探测器中心距离的增加而减少，并且在锥形光纤锥度比一定的条件下，改变锥形光纤长度，光强度随着长度的增加而下降，但通过锥形光纤的光子分布是一样的.这些研究结果对锥形光纤的设计和实际应用提供了重要依据.     

混沌实验数据处理及仿真

大学物理实验中混沌实验大多采用观察现象的方法进行,本实验采用蔡氏电路(Chua s circuit)产生混沌行为.在观察不同初始值条件下出现的倍周期分岔、阵发混沌、奇异吸引子等相图及现象的基础上,通过对采集数据进行处理,对负电阻伏安特性进行分段线性拟合,用功率频谱法、计算机仿真方法(龙格-库塔数值积分法)对混沌现象进行描绘,将实验数据与非线性方程组的数值解相结合,呈现出混沌现象的本质.     

强流脉冲离子束辐照靶及其喷发的数值研究

利用拟合实测的TEMP Ⅱ型加速器磁绝缘二极管(MID)电压波形及其焦点附近束流密度曲线，建立了Gaussian分布模型，据此计算了与靶作用的离子的能量及数量，采用Monte Carlo(MC)方法计算了沉积在靶内的能量.并以此作为热源，与流体动力学(HD)模型相结合，对不同的靶状态采用相应的状态方程，模拟计算了靶内压力演化情况; 同时对烧蚀产生的等离子体采用理想气体状态方程, 结合HD方程组, 模拟计算了喷发过程中压力的空间演化过程. 关键词： 强流脉冲离子束 Gaussian模型 HD方程 数值研究     

腔内倾斜扰动对正支共焦腔本征模式的影响

针对激光器出光过程中腔内像差扰动带来的光束质量下降问题，分析了正支共焦腔腔内不同位置倾斜扰动对本征模式的影响，采用数值迭代法计算了理想状态及腔内不同位置倾斜扰动状态下的本征模式、相位Zernike像差拟合系数.结果表明：对大小给定的扰动，相位相对敏感性从凸面镜到凹面镜逐渐增加，且腔模倾斜像差增大的同时还将导致离焦、像散等高阶像差的增大，因此作腔内补偿时，应使补偿平面尽可能靠近凹面镜.     

High-order accurate equations describing vibrations of thin bars

A method for deriving one-dimensional wave propagation equations in thin inhomogeneous anisotropic bars based on the mathematical homogenization theory for periodic media is used to obtain equations governing the longitudinal and transverse vibrations of a homogeneous circular bar. The equations are derived up to O(ε⁸) terms and take into account variable body forces and surface loads. Here, ε is the ratio of the bar’s typical thickness to the typical wavelength.     

不完全信息群体决策专家权重的集结

本文对于属性权重信息和属性效用信息都不完全的群体多属性决策问题，通过构造属性值区间和运用系统聚类分析法，对群体决策中的专家进行分类，并确定每位专家的权重．     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

聚丙烯孔洞铁电驻极体膜的电极化及其电荷动态特性

根据栅控恒压电晕充电组合反极性电晕补偿充电法的实验结果计算出铁电驻极体的极化强度.结果说明，伴随着薄膜内孔洞气体的Paschen击穿，该铁电体的极化强度随栅压增加而显著上升.利用上述充电方法和热刺激放电(TSD)谱的分析讨论了这类空间电荷型宏观电偶极子，及与其补偿的空间电荷热退极化的电荷动态特性；阐明了这两类俘获电荷的能阱分布，即构成宏观电偶极子的位于孔洞上下介质层内的等值异号空间电荷分别被俘获在深、浅两种能值陷阱内，而位于薄膜表面层的注入空间电荷则被俘获在中等能值陷阱中. 关键词： 反极性电晕补偿充电法 铁电驻极体 充电电流 热刺激放电